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ABSTRACT 



Parameter plane studies of an optimal second order 
regulator are presented. Emphasis is placed on the 
interpretation of the cost function and the sensitivity of 
the cost function to plant parameter incremental varia- 
tions . An analysis is made of cost function weighting 
factors and their effect on damping, speed of response, 
and cost. 
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I. 



INTRODUCTION 



The parameter plane study of an optimal regulator was 
undertaken with the purpose of finding an interpretation of 
the meaning of the cost function and of the weighting fac- 
tors as they relate to the observed physical behaviour of 
the regulator. 

Chapter one gives a brief description of the problem 
formulation; Chapter two describes some studies on the 
performance index; Chapter three deals with the study of 
the weighting factors, and Chapter four with some sensi- 
tivity aspects. Appendices A and B deal with the develop- 
ment of the general expression for the R matrix and the 
performance index. 

A. PROBLEM FORMULATION 

Suppose that initially the plant input or any of its 
derivatives is nonzero. Provide a plant input to bring the 
output or its derivatives to zero. In other words, the 
problem is to apply a control to take the plant from a 
nonzero state to the zero state. This problem may occur 

I 

where the plant is subjected to unwanted disturbances that 
perturb its output (e.g., a radar antenna control system 
with the antenna subject to wind gusts) . 

The desirable properties of a regulator system should 

be : 
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A. The regulator system should involve a linear 
control law. 

B. By definition an optimal system is one that mini- 
mizes a certain cost. 

To achieve property B, let us define a performance index 

oo 

PI = / (x'^QX + u'^Pu) dt 
o 

where Q is symmetric positive definite and R is a positive 

OO 

T 

definite matrix. / X QX dt come from the minimum integral- 

0 

square error problem. 

00 00 

/ E (X. ) ^ dt = / (x'^x) dt 

o i=l ^ o 

constituting a reasonable measure of the system transient 

OO 

T 

response, and / u Pu dt comes from the minimum energy 
o 

problem. 

The minimization problem, i.e., the task of finding 
an optimal control that minimizes the performance index, 
turns out to be achievable with a linear feedback law. This 
is the reason why the performance index includes a measure 
control energy. 

For the infinite-time regulator problem R(t) is the 
solution of the reduced Riccati equation 

R A + A*^R - R BP“^b'^R + Q = 0 
o o o o 

with an optimal control defined by 

u*(t) = -p"^(t)B’^(t)R^ x(t) 
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where 



A = 



0 1 
IP -P. 



, = 



~0 0 




“ 0 “ 


T 


" 




, B = 




, = 


0 K 


J- "P- 




_ K _ 




mm 



R = 
0 



^11 ^12 



r r 

12 22 



From Appendix A the general expression for the R matrix 



IS : 



^o = 



±p/q 



11 



/q 



K 



+ (±-jp-) (-P +/ P^+K^q22 ± 2Kq^^) 



:/q 



11 



K 






K 



p^ °22 

2 ^ 



ZR + 

K 



(l.D 



Uq(x) = -p“^b'^R^x 



if P ^ = 1, = [0 K] 



U^(x) = [±/q^i 



X, + 



2/q. 



'-P /v\ , --11 

T--^ JZ2-^ ^22 ^ 



K' 



K 1^2 



from which the system block diagram is: 
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Reducing to single block: 




“n " 



11 



( 1 . 2 ) 



C = 



2 . „2 



2v^ 



11 



2/K/qii 



p* + K‘(qj2 ± ) 



1 /_£^ • ^ ^22 



K/q^ K/q^ 



± 2 



( 1 . 3 ) 



B. OPTIMAL SECOND ORDER SYSTEM 




K and p are nominal values and and K2 are adjusted so 
that 
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I 



♦ 



(1.4 



= 



^ , K. = ^ 



'11 



2 K 



-p + /p- + K- 

\/ 0 Op 



2 *^22 /*^11 
± 2K ' 



11 



Oi/ P 



11 



and defined as in 1.1. 

The cost function for the optimal regulator is given by 
■^o = + 2r^jX^(o)Xj(o) + rjj x^ (o) 



‘11 



'11 



K J Pn 

o V ^11 



p2 * (522) + /5ll 

o O P^^ O '/ P^^ 



(0) 



2 nil 



X, (o) x_ (0) 



^ov/Pii 



K 



o >— 



■P + 
^o 



2 2 *^22 
p + K (-^) + 2K 



‘11 



o q 



11 



OU P 



11 



x2(0) 



And substituting the values of and gives 



'11 



^ (Po + KjK^) x2(0) 



2K 

IT 



^2 2 

+ x^CO) 



10 



II. THE PERFORMANCE INDEX 



A. GENERAL 

In Appendix B the general expression of tlie performance 
index for a second order system has been derived, giving, 
for initial conditions on X only 



(x^ (0) ] 



2 *^11 2 
(xf + (p + K K ) 

1 P-|, O 20 

(K" + — + K, K (K"+ -^) 

1 Pii K^Kq 1 o 2 p^^ 



.2 , ‘^ll, 



.2 , ^ 22 , 



2^2 <Po ^2 ^oM 2 (p- ■ 

J ^o 2 o 



( 2 . 1 ) 



where is the plant pole, 

and are the gains of the feedback paths, 
^11' ^^22' Pll weighting factors - 

Let = 1, q ^2 ~ 0/ ^ ~ 1600 and p^ = 30 

then ;r = 0.04002. 



[x^ (0) ] 



When q 22 = 0/ then equation (1.4) can be written as: 



K - ^ P * 
° ■'a ° 



which is a straight line on the K vs . p plane. Figure 2.1 
gives the values of K required to optimize the regulator. 

For example K = .02 is necessary to optimize a regulator 

. 2 

with a plant pole of 30 and = 2000. 
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Figure 2.1. Required values of K 2 to optimize the regulator. 
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B. STUDY OF J ON THE PARAMETER PLANE 

The cost function for the optimal system was given in 
section A by 



+ ( 2 . 2 ) 



1 [2 

substituting K 2 = ^t~P + Jv + 2k ] into 2.2 gives, for 
^^22 " ® 





+ 2K 



(2.3) 



Figure 2.2 shows the parameter plane for the optimal cost 
function obtained by varying K 2 when p changes , with K as a 
family parameter and = 1 . Lines of constant p and con- 
stant K have been drawn for clarity of presentation. 

The value of ^ corresponding to each point in the plane 
is obtained from 



2 

s 



+ 



(p + KK2 )s + KK^ 



0 



^ = 



p + KK 2 
2/KK, 



A plot of constant C curves on the , K 2 parameter 
plane for K = 16C0 and p = 30 is shown in Figure 2.3. 



From 2.1, for q 22 = 0 



KK 2 + (2p- 2KJ^)K2 + (K^ 



, 2 +"^ 11)2 _ 
Pll 



2pJ^ = 0, 



Choosing p = 30, K = 1600, = p^^^ = 1, the contour in 

Jl - K 2 plane is described by 



13 



i 



i 



♦ 




Figure 2.2. Loci of optimal cost function, = 1.0. 
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Figure 2.3. Constant C curves when K = 1600, p - 30 . 
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K 



-{.0375 - 2J^) 



± j '/(.0375-2J^)^ - 4 (.001953-. 0375J^) (2.4) 



2 



2 



and if K 



2 



constant = .02125 



128 




+ .56256 + K^j 



(2.5) 



Equations 2.4 and 2.5 are represented on the parameter plane 
in Figure 2.4. It can be readily seen that for = .02125 
the curve shows a minimum cost of . 04 . The same minimum 
cost is obtained keeping K 2 constant at a value of .02125 
for = 1. 
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Figure 2.4. Optimal cost function given or K 2 . 
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Ill': STUDY OF WEIGHTING FACTORS OF THE PERFORMANCE INDEX 

When it is desired to optimize the design of a control 
system, a cost functional must be established. This cost 
functional is an analytical expression, related to the 
system, which must be minimized by choosing parameters 
which will cause an extremal of the cost functional. 

Weighting factors used in performance indices deter- 
mine to a large extent the optimal system which results. 

A. SELECTION OF THE Q MATRIX 

Optimal control law determination by Athans and Falb [1] 
require that the weighting factors matrix be positive 
definite . 

Tuel [2] developed a canonical form for the weighting 
factor matrix Q which has the minimum number of parameters 
required in the performance index for the computation of the 
optimal control law. 

The selection of this Q matrix is a very important part 
of optimal design but there is very little guidance in the 
literature in the selection of the weighting factors. 

In any regulator the output or "controlled variable" is 
of primary interest and certainly must be weighted, being 
never negative. Then 

‘3ll " 0 

The velocity (X 2 = need not be weighted at all 

but certainly any weight attributed to the velocity will 
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alter the dynamic response and may give desirable features. 
Thus • 



But we should normally expect q ^2 

B. DERIVATION OF EXPRESSION FOR STUDY IN THE PARAMETER PLANE 
For the second order system the characteristic equation 



IS 



s + (p + K..K )s + K K, = 0 
'^o 2 o o 1 



where 



Let 



then 



1 V p 



11 



11 



K = - 1° + -L 
^2 K K 

o o 



/ 2 2 *^22 t 

'' o o p^j^ oV ] 



11 

’ll 



then 



p + K K- = 
^o o 2 



2 . „2 ^22 



p + K 
^o op 



11 



± 2K = 2?W 



o V p 



11 



n 



K K, = CO = K 



o 1 n o ^ p 



^11 



11 



K 



■oVpii 



= y 



^22 

’ll 



= X 



2 2 2 2 
4C 0 ) = p + K X + 2y. 

n o ^ 
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It can be noted that the constant C curves are straight 

lines , but the location of the family for 0 < C < 1 depends 

on K and p , i.e., the effect of q, , and q__ can be shown 
o o 11 ^22 

only for a specific plant. 

Figure 3.1 shows the parameter plane for the plant with 

K = 1600, p = 30. 
o *^o 

From section B the cost function for the optimal system 

was : 






± 2K 

Pll 




( 0 ) 



+ 




(0) x^ (0) 



Let 





+ K 

o 



+ 2K 

Pll 



0 




x^(0) 



J 

O 



Pll [x^ (0) ] 



■'oh 



K 

o 



Pll 




+ K 




) ± 2K 

o 



Pll 



= yVp^ + K^x ± 2y 
-^*^0 0 ^ 

For the optimal system: 

= .04, K = 1600, p = 30. 
1 o ^o 
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Figure 3.1. Constant curves on - q^2 parameter plane 



vs . 



A constant cost curve is shown in Figure 3.2 in the 
2 

0 )^ parameter plane. 




C. EFFECT ON DAMPING AND COST, STEADY STATE ACCURACY 
AND SPEED OF RESPONSE 

In order to study the effect of the weighting factors 
on damping in the parameter plane it is necessary to use 
1 . 2 and 1.3. 



0) 



2 

n 




= z 



= iV + 



K 



2 

^11 



± 2 



let 



then 



^22 = ”'3ll 






+ Wz ± 2 
z 



(3.1) 



From this equation the following expression can be obtained: 
p^ = -z(±2 - 4?^) -Wz^ (3.2) 

As K is defined in Section A as: 



K.. = 



.. E + 1 

K K, 



p + K 



2 

Pll 



± 2K, 



'll 

’ll 



where 



K 



1 



+ 
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Figure 3.2. Cost function curve when K^=1600, p^=30. 
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Then, for q = 0 the choice of + sign will require nega- 
tive feedback, while the - sign will require positive feed- 
back . 

For negative feedback 3.2 can be rearranged to: 



Then, as z must be positive, C must be greater or equal 
to .707 noting that if C = .707 then p = 0. 

In Figure 3.3 the parameter plane for W = 0 has been 
drawn for different values of C and in Figure 3.4 for 
VI ^ 0 . In both figures, curves of constant cost function 
have been superimposed. 

It can be noted: 

1. Minimum C for optimal regulator (second order) with 
negative feedback is .707, and this occurs only when W = 0. 

2. For W > 0, C > .707. 

3. As W increases, the system rapidly becomes over- 
damped . 

4. If p is large the system will be very heavily 
damped, no matter what the gain is. 

5. For W ^ 0 if the gain is raised, the system C goes 
thru a minimum and then increases with gain. 

6. The larger W, the more velocity feedback is required. 

For plotting of cost functions , K has been kept constant 



X 



for W = 0. 



2 - 4C 



at 1600 




stant at 1 and K varied in Figure 3.3b. 
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From these two figures the following results can be 
established: 

1. For K constant, the cost function increases with, 
increasing weighting factor. 

2. For constant weighting factor the cost function 
decreases with increasing gain. 

Figures 3.4a thru 3.4d show the parameter plane for 

W 0 and it can be readily seen that for any p and ^ two 

values of can be obtained, both having different cost. 

/qii 

In Figure 3.4a and 3.4b, >/ has been kept constant 

Pll 

at 1 but for the cost analysis a value of W = 1 has been 
used in 3.4a and W = . 5 in Figure 3.4b. A comparison of 
these curves shows a larger cost for larger W, and the cost 
decreasing with an increase in gain. 

On Figures 3.4c and 3.4d the gain has been kept constant 
at 10 and the weighting factor varied. These figures 

show a decrease in cost for a decrease in and a slight- 

ly smaller cost for larger W. The value of p^^^ has been 
kept constant for these and the following cases. 

When using positive feedback any value of damping can 
be used and the results are shown in Figures 3.5 thru 3.7d. 

In Figure 3.5, K has been kept constant at 1600 and 
varied and in Figure 3.6, q^^ has been kept constant and 
the gain varied. 

The loci of equal cost function follows different 
patterns in both figures but in general an increase in cost 
follows a decrease in gain. 
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Figure 3.4. Constant ^ curves for different values of 
W, negative tachometer feedback. 
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Figure 3.4b. Constant cost curves for W = .5 
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1, negative tachometer feedback. 



J=1.0 



1.25 

.5 



I P 
2.5 

and 
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Figure 3.4c. Constant J curves for W = 1.0 and K = 10, 
negative tachometer feedback. 
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Figure 3.4d. Constant J curves for W = .1 and K = 10 , 
negative tachometer feedback. 
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'igure 3.5. Constant c; curves for W = 0 and positive 
feedback for constant J curves, K = 1600. 
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A comparison between Figure 3.3 and Figure 3.5 shows 
that -for negative feedback values of ? ^ .707 are restricted 
to the first quadrant while for positive feedback any value 
of C is permissible. 

Figures 3.7 show the parameter plane for two different 
values of W and for different C* 

A comparison between negative and positive feedback 
shows same results as for W = 0, i.e., any value of C is 
permissible in the first quadrant for positive feedback. 

Figure 3.7a and 3.7b compare the effect on cost of 
different values of W, keeping q^^^ constant and varying the 
gain which results in larger cost, which for p = 1.5 is 
0.75 for W = 1 and 1.25 for W = 2. 

Figures 3.7c and 3 . 7d compare effect on cost of dif- 
ferent values of W, keeping the gain constant at 10 and 
varying For a small value of p the cost is smaller 

for W = 1 than for W = 2 and for any there is a minimum 

cost, which for w = .7 is about 0.19 for W = 1 and 0.3 

n 

for W = 2 . 

Keeping C = 1.0 and K constant at 10, a curve of cost 
against frequency has been drawn and the results shown in 
Figure 3.8. It can be seen that the cost follows the shape 
of a parabola increasing very fast as the frequency 
approaches 2 . 

In order to study the effect of the weighting factors 
on the steady state accuracy it is convenient to refer to 
the equivalent block diagram of Chapter I. It can be seen 
that 
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Figure 3.7. Constant J curves for W ^ 0, positive 
tachometer feedback. 



36 




±r_ I I I L 

1 2 3 4 5 



Figure 3.7a. Constant J curves for W=1 and 
positive tachometer feedback. 




= 1 , 
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Figure 3.7b. Constant J curves for W=2 and 
positive tachometer feedback. 




= 1 , 
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Figure 3.7c. Constant J curves for W = 1 and K = 10 , 
positive tachometer feedback. 
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Figure 3.7d. Constant J curves for W=2 and K=10, 
positive tachometer feedback. 
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Figure 3.8. parameter plane for C=l-0 and K=10. 
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which will depend only on and the parameter plane is 

shown in Figure 3.9. For small values of the steady 

state error is also very s mal l , increasing rapidly as q^^^ 



ing 1.0 in the infinite. 

For the speed of response study it is convenient to use 
the definition of settling time, i.e.. 



Several cases should be considered: 

1. Negative feedback for the velocity. Curves of 
constant settling time have been drawn in the parameter 
plane for two different weighting functions, W=0 and W=2, 
the results being shovm on Figure 3.9a. It can be seen that 
the settling time decreases with increasing W. 

2. Positive feedback for the velocity. Constant 
settling time curves for the same parameters as in the pre- 
vious case are shown on Figure 3.9b. This curves show that 
the settling time decreases with increasing W but faster 
than in the negative feedback case. 

3. Negative feedback for the position. The results 
are as in Case 1. 



increases to a value of 




and then slowly approach- 




n 
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Figure 3.9. 



Steady state error vs. 




parameter plane. 
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Figure 3.9a. Constant settling time curves for negative 
tachometer feedback. 
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Figure 3.9b. Constant settling time curves for posi- 
tive tachometer feedback. 
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4. Positive feedback for the position. Using Routh 
criteria of stability it can be seen that this case gives 
a pole in the right hand plane, meaning that the system will 
be unstable and the regulator will not regulate at all. 



46 



IV. SENSITIVITY IN THE OPTIMAL REGULATOR 



A. GENERAL 

An optimal design guarantees minimum cost for the regu- 
lator from any initial conditions. 

If any parameter of the regulator deviates from the 
value used in (or determined by) the optimal design, then 
the cost of regulating is increased, i.e., the performance 
of the system is not optimal. 

Sensitivity is a word used to describe the rate at 
which some characteristics of the system deviates from a 
reference value as a function of a parameter change. Thus 
various types of sensitivity could be defined: 

a. Root sensitivity 

b. Bandwidth sensitivity 

c. Steady state accuracy sensitivity 

d. Settling time sensitivity 

e. Rise time sensitivity 

f. State sensitivity 

g. Cost sensitivity 

State sensitivity is a measure of the deviation of a 
dynamic state from the values it would assume if the system 
was optimal. 

Cost sensitivity is some measure of the deviation of 
the cost from the optimal cost. 

Ultimately all of these various sensitivities relate to 
the same basic characteristics of a linear regulator, i.e.. 



they indicate changes in root location (possibly excepting 
steady state accuracy sensitivity) . However there must be 
some basic differences in these various sensitivities in 
the sense that some are vector quantities and others are 
scalars, i.e., root, bandwidth, settling time and rise time 
sensitivities all indicate pretty clearly the direction in 
which a root has moved, while cost sensitivity only indi- 
cates the magnitude of the root motion. 

From a design point of view there may be some advan- 
tages to the cost sensitivity, i.e., one can usually accept 
dominant root location within a specified s-plane area, 
and a defined area on the parameter plane normally maps 
into a dominant area on the s-plane. Correlation with 
other performance criteria may very well be required. 

Cost may be evaluated at any point on the parameter 
plane and thus constant cost contours can be obtained. 

Cost sensitivity in a macroscopic sense is then just the 
difference between two costs divided by the increment in 
parameter value between them. 

Cost sensitivity can also be evaluated at a point and 
normally the point of greatest interest on the parameter 
plane would be the point at which the parameter assume 
optimal values. Change in cost from this optimal value 
(for a small change) can be computed by: 

a (paraLter) ^parameter = AJ 

J^+AJ=J^ 
opt at new point 
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This however depends on initial conditions values, 
and the permissible range of A parameter depends on the 
system and on the cost function. 



B. SENSITIVITY OF THE COST FUNCTION AND SENSITIVITY 

OF THE OPTIMAL COST FUNCTION 

In the research, the curves of constant J on the parame- 
ter plane are obtained using an expression for J which is 
derived in terms of p, K, , K2 , where p and K are plant 
parameters and and K2 are values for the state feed- 
back loops . In order to use the expressions the system is 
optimized at some chosen p and K, and values of and K2 
are computed for the system thus optimized. When J is set 
to a non-optimal value and K2 remain fixed at the opti- 
mal values and p and K are computed. The system is not 
optimal for the new p and K. 

Sensitivity of the optimal system could be defined by 
taking the derivative of the cost function, i.e.. 



S 

op 



9 J 



8 p 



p=p , K=K 
^o o 



S 



J 

oK 




K=K 

o 



However the equation for J optimal can be written in 
at least two forms : 

A. The form is used to compute J curves, which con- 
tains the feedback gains and K2 as constants. 

B. A form obtained directly from the R matrix such that 
and do not appear as symbols, the cost being 
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expressed entirely in terms of p, K, <3^1' ^ 22 ' ^11* 

Both forms A and B should give the same numerical 

answer when p = p , and K = K . . However their deriva- 
^ '^opt opt 

tives at the optimal point should be different, because 
form A describes a system that is not optimal when p and 
K are changed from their optimal values, while form B 
presumably expresses the condition whereby changes in p 
and K automatically result in changes in and so 
that the system remains optimal at the new p, K values, 
but of course with different value of J. 

Under these conditions the choice of a definition for 
sensitivity becomes a matter of concern: 

Should sensitivity be defined on the basis of jj] 
from the value it has when p and K are at their specified 
optimal values? or should be defined on the basis of the 
deviation of jjj at point p, K, from point p,K? 

By this last statement it is meant: 

If and K 2 are evaluated at p, K, then 

= f<Po' >'l' >'2> 

At some new set of plant parameters values p^, K , 

3. a. 

A = *<Pa' ''a' >"1' ■'2> 

which corresponds to the curves drawn on the parameter 
plane, but the system is optimal at p , K . If the system 

d. d. 

were optimal at p , K , K, , K- , would have to be deter- 

3 3 _LS Z3 

mined and the cost function would become 

■^la = *<Pa' "a- ^la' '^2a' °Pti">al- 
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C. SENSITIVITY OF COST FUNCTION TO PLANT PARAMETERS 
VARIATIONS 

If is adjusted for the optimal value, then 



= \/ - 
1 V p 



11 



11 



and 



^11 



.2 . ^22, 



2kI + (p + KjK)^ + K^K(K‘ + 



-2K^K2(p + K^K) 
2K, 



x^(0) 



2(p + K 2 K) 



— Xjl(O) X2(0) 



2Ki 2 *322 

+ KI + 



K 



2 p 



11 



4io) 



2(p + K 2 K) 



If K 2 is also adjusted for the optimal value then 



2K p^ 2K^ 

— Kl + — ^ 
p-.-, 2 K, 



^22 „2 . *^0 ""'1 



'11 



K 



and 



'11 



2 2K 

2 Kt + — - 

1 K 



(p + K 2 K) + K^K(2K2 + 



2K p 2K. 



K 



K 



-2Kj^K2(p + K 2 K) 
2K, 



(0) 



2(p + K 2 K) 



+ 



— Xj^(0)Xj(0) 



2K, „ 2K„p 2K, 



x2(0) 



2 (p + K K) 
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This is the cost function assuming that and K 2 are 
adjusted to optimum for the nominal values of system 
parameters, p^ and . This result may be illustrated more 
graphically by rearranging as follows: 



Pll 



2r2 ( 1 - + 2K K (j^ p - p) 

o o 



2K 

-f (p + K^K)^ 

2K 

+ (0) X 2 (0) 



x^ (0) 



2(p + K 2 K) 



+ 



2Ki (i - i ) + 

o 



2K_(p^ + K.K ) 
2 O 2 o 



K, 



x^(0) 



2(p + K 2 K) 



The sensitivity becomes: 



J - J 



921 = 



'11 



2 K ^Po " P^o 

K (1 - f) + K K (^5^ £) 

o o 

p + K 2 K 



K K 

+ ^ (p + KjK) - X <Po + 



X^ (0) 



+ 2 



K K 



(0) X 2 (0) 



K (^ - 1 -) 

1 'k K ^ 

— + — 

p + K«K K 



p + K„K 
^O 2 o 

p + K^K 



- 1 



x^ (0) 



(4.1) 



'^1 



J - J 



921 =. 



11 



K K 

^ (K - K) +4-2 (Kp - pK ) 

-5 2 + K, (| - 4) 

p + K 2 K ^ ^ 



x^(0) 



52 



+ 2K, 



1 _1_ 
K " K 



(0) X 2 (0) 



+ - 



Kq - K 

^ - P) + K2(K^-K) 



P + K2K 



k^(p+K2K) 



XjCO) 



Note these become zero when K K , p p 

o ^ ^o 



J - J 



opt _ 



ll 



k; 



+ 



(p+KjK) 

2K, 

<V") Xl«>>=<2«» 



K K K 

'"o-"* ^ IvMpo - ks; J 



x^(0) 



+ 



Ik, k'z ] 




12 

KK K 


^ ir <Po-p> 


0 

0 


0 



x^ (0) 

(P+K 2 K) 



For small changes where 



K = K + 6K 
o 



K K 



p = p +6 



P = P, 



6J 



11 



K 



+ 



10 o 2 o 
1 



(1-6K) 



K^K2 



K. 



K 



Ip +K„K K 
o ^^o 2 o o 

2K, 



(6K)p^ - K^(6p) 



K 



— (6K) (0) X 2 (0) 



x^ (0) 





' 1 


(k^ kA 


K„ 


+ . 




1 2 ' 

,^2 K 


(-6K) + ^ (-6p) 




1 


1^0 


1 0 



x2(0) 

Po+KjK 



6J 


-6K 


r 2 2 

p K^+KfK 
^0 1 1 0 


^lPo-^^1^21 


-6p 


0 

M 

1 




Pll 1 


k2 
^ 0 


P +k„k 
^0 2 0 


K 

0 


_"o<Po+W. 





x^ (0) 
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2K 



1 




(6K)Xj^(0)x2(0) 



\ 



+ 




(6p)> (0) 
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From 4.1 if p = p^ = 30 and K deviates from the optimal 
value of 1600 a curve can be represented in the parameter 
plane and is shown in Figure 4.1. Also shown is the curve 
for K = = 1600 and p deviating from the optimal value. 
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% 





Figure 4,1, Sensitivity of cost function to plant 
parameter variation. 
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V. CONCLUSIONS 



RESULTS 

An optimal regulator can always be obtained using the 
parameter plane method for any given plant. Both negative 
and positive feedback can be used to achieve the desired 
results but only negative position feedback should be used 
to obtain a stable system. Positive position feedback 
gives an unstable system with no regulation at all. 

For a given plant an optimal cost function can be 
obtained and the values of the feedback path gain calcu- 
lated using the curves of Chapter II. 

When a given frequency and damping is desired, the plant 
parameters needed to achieve the results can be obtained 
for different weighting functions with the use of figures 
of Chapter III, or for a fixed plant the performance and 
cost be obtained by the use of the same figures. 

For plant parameters variations the sensitivity of the 
cost can be obtained using the figure in Chapter IV. 
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APPENDIX A 



GENERAL EXPRESSION FOR THE R MATRIX 



From the reduced Ricatti equation: 

A*^R + R A - R BP'^b'^R + Q = 0 
o o o o 



where 



A = 



0 1 
0 -p 





0 


0 




"n 


= 






B = 






_1 


"P_ 


/ 


_K_ 



(A.l) 



b'^ = [0 K] 



R = 
o 



^11 ^12 



T r 

12 22 



Substituting into A.l 
I 0 

-P 



^11 ^12 


+ 


^11 ^12 


Jl2 ^2^ 




Jl2 ^2^ 


^11 ^12 


+ 


qil 0 


_fl2 ^ 22 ^ 




1 ^2^ 



0 1 
0 -p 

= 0 



^11 ^12 



T T 

12 22 



0 K 















0 


0 


+ 


0 


^11-P^12 




^ ^12 ^ ^12^22 
t ,2 ,,2 


3i~^^22 


^12~P^22_ 




0 


^12"P^22_ 




K r^2 K r22 



^11 
0 q 



0 

22 



= 0 



v2 2 ^ 

-K r^2 + ^11 



= 0 



^ll"P^12 ^ ^12^22 ® 



(A. 2) 
(A. 3) 
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0 



(A. 4) 



rii-pri2-K r^2^22 



2 2 

2 (ri2-pr22)-K = 0 



22 ^22 



(A. 5) 



From A. 2 



^11 






12 



K 



r = + 

12 - K 



11 



Substituting in A. 5 



K 



,2 2 



2pr22 K T22 + <^22 ® 



Since Q must be positive definite must be positive 

but both signs will be used with the radical. 

. 



+ 2pr,, - q 



22 



22 



= 0 



K 



22 



-2p ± sj 4p^+4K^ (q 22 ± 

2k2 



2 /g 



11, 



K 





+ 4 r 2 (^22 ± 



2r2 



2 ,/qii 

K 



As r 22 niust be positive only the positive sign outside 
the radical can be used, and the negative sign inside the 
radical only if 

2 2 J 

'922 - — > i ° 



i.e. for 





> 


2 '^ll 


p22 




K 



Then, both signs will be carried with the indicated restric- 
tion implied. 
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From A. 3 and A. 4 





K 



y?^+K^q22 



± 2K/5^ 



Thus the general expression for the R matrix is; 



/q 



11 



K 



± 2K/q 



11 



- 1 - 



/q 



11 



K 



R 



o 




11 
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APPENDIX B 



GENERAL COST FUNCTION FOR A SECOND ORDER SYSTEM 

The following derivation was originally worked by Dr. 
Sidney R. Parker and has been included because its extensive 
use in this thesis and not being published in any paper. 

It has been shown that for the general second order 
regulator 



J = / (x’^Qx + u'^pU) dt 



where 



Q = 



qil 0 



'22 



/ P “ P]^]^ snd U = 



K 1 K 2 



X, 



If 



T T T 

X Qx + U pU = X MX 



(B.l) 



Then 



T T 

X Qx + U pU = 



Xi X 2 



Xi X2 



qil 0 



^22 






X, 



Xi X 2 




Ki K 2 



X, 



^11 









\ 




0 

^ 22 _ 


4- 


Ki K^K2 

V 2 ''2_ 


Pll) 


^1 

^ 2 _ 






qil+KiPjl 

Pll^l^2 



P11V2 

^ 22 ‘’'^ 2 Pll 



X, 
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Then from B.l 



mil 


= ‘3ii 


+ 


^iPll 


"’12 


"’21 


= 


Pll^l^2 


m 2 2 


^^22 


+ 


^2^11 



Then : 



T 

X MX 



= x^] 



"’ll "’l2 

"’12 "^22 



If: 



then 



- m ^ 2^1 ^’^ 12 ^ 1^2 "^ 22^2 



-Xit -^2t 

^1 ~ ^‘l^ “2 ^ 



~^1^ 

Xj = -X^o^e - XjOjS 



and 



and 



then 



2 r 2 ~^^1^ ~ ^'^l 

1 ~ "’ll[^‘^l® ''’ ^a^a^e 



= 



/ m^^x^dt = 
o 



-(A^+X2)t 2 “2^2^* 

+ a2e 



^°‘l‘^2 “2 



2 X 1 ^ 1+^2 ^^2 



(B.2) 



"’ 22^2 "^22 



. -2X t -(X +X )t » -2X t 

^1*^1^ + 2X^X2a^a20 +X2a20 



^ "’ 22 ^ 2 ^^ "^22 

o 



""^l^^l , ^^1^2°‘l°‘2 , ^2^2 

-j- . ^ “ ”1" 



2X- 



X 1 +X 2 



2X. 



"’ 12 ^ 1^2 = "2’^12 



then 



2 -2X,t 2 “^^ 2 ^ 

Xitt^e + ^2*^2^ '*' “2^062 ^^l''’^2^ ^ 



(B.3) 



-(X^+X2)t 



/ ^"’12 

o 



2 A 2 -, 

^l“l , . ^22 

+ a,a„ + 



2X. 



'1^2 2X. 



(B.4) 



61 



And 



J = B.2 + B.3 + B.4 



= m 



11 



a 2 
1 

.2Ai 



"T " + 



X1+A2 



2A 



2J 



- 2m 



12 



'^ 1*^1 ^ ^ 1^2 ^^ l '*'^ 2 ^ ^ ^ 2°^2 



2X. 



A 1 +A 2 



2A 



2 J 



+ m 



22 



,22 ^ 2 , 2 

*~^l‘^l ^ ^^1^2'^l‘^2 ^ ^ 2 ^ 2 "’ 



2A. 



A 1 +A 2 



2A 



2 J 



To find and A 2 let: 



(B.5) 




X-, 



^ + (a+K2G)s + K^G 



w = /K,G 
n 1 
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Then 



= a + K 2 G 



? = 



a + K 2 G 



2 



Let 



2 2 
s Xt + 2^00 sx, + 00 X, = GR 
1 n 1 n 1 



Xi = Xi , X2 = x^ 



x^ = -2^00 x^ - (JL) Xt + GR 
2 n 2 n 1 



• 




Xt 




1 


_ 


• 




_^2_ 





-(I) -2^w 

n n 



X. 



X , 



when R = 0 



X 



= [sI-A 



-1 



x(0) 

'V 



sl-^ = 



-1 



0) s + 2Cw 

n n 



[si-a] 



-1 



2 2 
S +25(ji) S + (j0 

n n 



s+2Cw^ 1 



-o)n 



^1 


1 


S+ 2 


1 




x^ ( 0 ) 


^2_ 


”2 2 
s +2CU) s+oj 
^ n n 


2 

-Cl) 

n 


3 




3‘°L 



X. = 



(s+2Cw^)x^ (0)+X2 (0) 

2 2 
s +2?w s+o) 
n n 



(s+Xj^+A 2) Xj^ (0)+X2 (0) 

(s+X^) (s->X2) 



^2 



-oj^Xi(O) + 3 X 2 ( 0 ) 

2~Z 2 

S +2^0) S->03 

n n 
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where 



“l , “2 

X, = — -r + 



1 s + Xj ^ ^"^^2 



Ai = cu^ + 

^2 ' '“n - “n 



a 



(s+2C(i)^s) (0) + (0) 



s+ 



s=- A. 



(-A + 2 ^ 0 ) )x, (0)+x^(0) 
1 n 1 2 

^2~^l 



a 



X^Xj^(O) + x^iO) 
^2~^1 



Similarly : 



“2 = 



X^Xj^(O) + X2(0) 

- (A2“A^) 



Note that ; 



(«! + a^) = x^(0) 



From B . 5 : 



'11 



= A 



r- 2 2 

a, a 

■^ + 2 V - 



2A- 



A 1 +A 2 2 A 2 



^1^2 j“i+2a^a 



2'^“2 



+ B 



r 2 _ 

a, 2a,a„ 

^ ^ ^ 



a. 



2 -I 



L^A]_ ^1 2 ^2 J 



where 



A = Kr + 



2 . ^11 



1 P 



11 



B = k:: + 



2 . ^^22 



2 P 



11 



(B.6) 



(B.7) 



(B.8) 
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Substituting B.6 and B.7 into B.8 and rearranging 



J 


f 

A 


^2 2A^A^ A^I 


Pll I 


1u2-Xi)2 


2X^ ^^2 



B 



^2 






2 “1' L^“l 



2X 



-1 



^^ 1^2 '^1 



2 



2X 



- ^l'"2 



(0) 



‘ (^2 ^1^ 



■^2 

^1 ^2 



-1 

2 J 



+ 



B 



a^-\) 






2 (X^+X^) 



-1 

‘2 J 



(0) X 2 (0) 



A 


- 

1 _ 2 1 




2X^ (Aj^+A2) 2 A 2 



B 



(A^-X, ) 



2X 



-1 



2 '^1' L"-'^l 

Substituting 

A- = ^(0 + u Vc^-1 = a + 

1 n n V 



x-i+x-i 



2X 



-1 
2 J 



x^CO) (B.9) 



X- = cw + w v/c^-1 = a - 

2 ^ n nV 



\ -i 2,2 

^1^2 ~ ^ 



into B.9. 



'11 



4b' 



(a-b)^ 
2 (a+b) 



2(a^-b^)_^ (a+b)^1 



2a 



2 (a-b) 



J 



B 



4b' 



-K 1 K 2 



(a-b) ^ (a+b) 2 (a^-b^) (a+b)^(a-b) 

2 “ 2a 2 



(0) 
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4b' 



a-b _2 ^ a+b 



a+b 



a-b 



(0) X 2 (0) 



A 



4b' 

B 



2 + 1 



2 (a+b) 2a 2 (a-b) 



a+b _ 2 (a^-b^) ^ a^b 



4b^ L ' 

Adding and rearranging 



2a 



X 2 (0) 



'11 



A 

4 



t- 2 , 2 , 2 , 2. 

5a -b I B(a -b ) 

, 2 ^2> i ^ 4a ^1^2 

i (a -b ) J 



( 0 ) 



+ 



2 ,2 

a -b 



x^ (0) X 2 (0) 



h + ^ 

4a(a2-b2) 



x^(0) 



Substituting 

a = Cw 



n 



= “r.y^ 



into B.ll 



Pll 



A(1+4C^) + Bu)2 - 4^0) K^K 



n 



n 12 4cw 



xj(0) 



n 



+ -y X. (0) X 2 (0) 
n 



x^ (0) 



^[4 



But 



A = K, + 



‘11 



1 Pll 

B = k 2 . m 

2 Pll 



w = K,K 
n 1 



(B.IO) 



(B.ll) 



66 



Then ' 



= p+K^K 



'11 



2 ^11 
Kf + 

1 Pll 



^11 

(K^+ (p+K K) ^ 

+ Itk 



KI+ 



2 P 



22 

111 



- 2K^K2(p+K2K) 



( 0 ) 

2(p+K2K) 



2 *^11 
k:+ 

1 Pll 

+ — — x^(0)X2(0) 



2 '^ll 
K,+ 

1 Pi 



'11 



K^K 



+ k::+ 



2. ^22 



2 P 



11 



x^(0) 

2 (P+K 2 K) 



(B.12) 
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